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We present the details of a method for conducting a targeted, coherent search for compact binary 
coalescences. The search is tailored to be used as a foUowup to electromagnetic transients such as 
Gamma Ray Bursts. We derive the coherent search statistic for Gaussian detector noise and discuss 
the benefits of a coherent, multi-detector search over coincidence methods. To mitigate the effects 
of non-stationary data, we introduce a number of signal consistency tests, including the null SNR, 
amplitude consistency and several tests. We demonstrate the search performance on Gaussian 
noise and on data from LIGO's fourth science run and verify that the signal consistency tests are 
capable of removing the majority of noise transients and the search gives an efhciency comparable 
to that achieved in Gaussian noise. 



I. INTRODUCTION 

There has been excellent progress towards gravita- 
tional wave astronomy over recent years. The first gen- 
eration of large scale gravitational wave interferometers 
reached unprecedented sensitivities and have undertaken 
extended science runs. The U.S. Laser Interferometer 
Gravitational- wave Observatory (LIGO) .Ij, the French- 
Italian Virgo [g and the German-British bEO600 [3| de- 
tectors now form a collaborative network of interferom- 
eters. The data from these detectors has been analyzed 
for gravitational waves from compact binary coalescence 
[1], stochastic background [S], unmodelled burst and 
pulsar [7] sources. LIGO's sixth science run (S6) and 
Virgo's second and third science runs (VSR2 and VSR3) 
ended in October 2010 and yielded the most sensitive 
data yet taken; the analysis of this data is ongoing. In 
the meantime, the detectors are being upgraded to their 
advanced configurations [8l410j. with the expectation of 
a ten fold improvement in sensitivity. With these sen- 
sitivities, it is expected that gravitational waves will be 
observed regularly |llj . Furthermore, with a proposed 
advanced detector in Japan [12], a possible detector in 
Australia [T3], and 3rd generation detectors on the hori- 
zon [Uj, future prospects are promising. 

As the gravitational wave community matures it is es- 
sential that a relationship is built between gravitational 
wave (GW) and electromagnetic (EM) astronomers. The 
GW emission from a source is likely to provide comple- 
mentary information to emission in various EM bands, 
and a joint observation is significantly more likely to an- 
swer outstanding astrophysical questions. Already this 
relationship is beginning to mature. A number of EM 
transients have already been followed up in GW data 
[T5HT7] . Additionally, infrastructure is also being put in 
place to allow for EM follow-up of G W observations [TH] ■ 

Compact binary coalescences (CBC) are one of the 
most promising sources of gravitational waves, and also 
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an ideal candidate for joint GW-EM astronomy. During 
the late stages of inspiral and merger, a compact binary 
emits a distinctive, "chirping" gravitational wave signal. 
Furthermore, CBCs containing at least one neutron star 
(NS) are expected to emit electromagnetically. Specifi- 
cally, binary neutron stars (BNS) and neutron star-black 
hole binaries (NSBH) mergers are the preferred progeni- 
tor model for the short gamma-ray burst (GRB) [roil20j. 
It is also possible that these mergers will be observable 
electromagnetically as orphan afterglows [19] , optical [2V 
or radio transients [22] • Since GRBs are well localized 
both in time and on the sky by EM observations, the cor- 
responding GW search can be simplified by reducing the 
volume of parameter space relative to an all-sky, all-time 
search. Targeted searches for CBC waveforms associated 
to short GRBs were performed using data from LIGO's 
fifth science run (S5) and Virgo's first science run (VSRl) 

[TgdS]. 

In this paper, we introduce a targeted, coherent search 
algorithm for detecting GW from CBC. This targeted 
search is designed as a follow-up to EM transients, and in 
particular GRBs. Previous searches for this source have 
made use of a coincidence requirement — namely that a 
signal with consistent parameters is observed in two or 
more detectors in the network. The analysis introduced 
here makes use of the data from all operational detec- 
tors, and combines the data in a coherent manner before 
matched filtering against CBC template waveforms. In 
a coherent analysis, it is straightforward to restrict the 
signal model to only two independent polarizations. This 
allows for the rejection of incoherent background noise, 
and consequently increases the sensitivity of the search 
if more than two detectors are operating. 

The data output by gravitational wave interferometers 
is neither stationary nor Gaussian, but is contaminated 
by noise transients of instrumental and environmental 
origin. This makes the task of doing data analysis a com- 
plex one, and matched filtering alone is not sufficient to 
distinguish signal from noise. Regardless of whether a co- 
incident or coherent search is performed, the most signif- 
icant events by signal-to-noise ratio (SNR) would always 
be dominated by non-Gaussian transients, or "glitches" , 
in the data. A significant effort goes into understand- 
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ing the cause of these ghtches [53] and removing times of 
poor data quahty from the analysis. While these efforts 
greatly reduce the number of glitches they cannot remove 
them entirely. Therefore the analysis must also employ 
methods to distinguish signal from noise transients. In 
previous CBC searches, signal consistency tests [24l ES] 
have proved very effective at removing the non-Gaussian 
background. We extend these tests to the coherent anal- 
ysis described in this paper and demonstrate their con- 
tinued effectiveness. In addition, coherent analyses natu- 
rally lend themselves to multi-detector consistency tests, 
such as the null stream [26 . We describe a number of 
such consistency tests for this templated CBC search and 
again demonstrate their efficacy. Finally, we are able to 
show that the various signal consistency tests are sufR- 
cient to remove the majority of non-Gaussian transients 
and render the search almost as sensitive as if the data 
were Gaussian and stationary. 

The layout of this paper is as follows. In section [llj 
we describe the formulation of a targeted coherent trig- 
gered search for CBC signals. In section [Ill| we discuss an 
implementation of the null stream formalism and other 
multi-detector consistency tests. In section IV we de- 
scribe a number of tests than can be applied in a 
coherent search to try to separate and veto glitches. Fi- 
nally, in[V]we outline an implementation of a targetted, 
coherent search for CBC and present results on both sim- 
ulated, Gaussian data and real detector data taken from 
LIGO's fourth science run (S4). 



II. COHERENT MATCHED FILTERING 



In this section, we describe the coherent matched- 
filtering search for a gravitational wave signal from a co- 
alescing binary in data from a network of detectors. We 
restrict attention to binaries where the component spin 
can be neglected. The description is primarily tailored 
towards searches where the sky location of the gravita- 
tional wave event is known a priori, as is appropriate 
when performing a followup of an EM transient such as 
a GRB jlH [19] . Finally, since all previously published 
CBC search results [H [17| [55] have used a coincidence 
search between detectors, we compare the coherent anal- 
ysis with the multi-detector coincident analysis. 

The coherent analysis for coalescing binary systems has 
been derived previously using a similar method in |291 - 
[55] . Our presentation makes use of the J'-statistic for- 
malism, introduced in [33j . This was originally defined 
as a method for performing searches for continuous wave 
searches and has been regularly used for this task (see for 
example [34]). It was noted in that the J-"-statistic 
and the multiple detector inspiral statistic derived in [30] 
are similar and the J^-statistic was adapted to searches 
for CBC signals in [36j . 



A. The binary coalescence waveform 

The generic binary coalescence waveform depends 
upon as many as seventeen parameters. However, we 
restrict attention to binaries on circular orbits with non- 
spinning components. This reduces the parameter space 
to nine dimensions: the two component masses Mi and 
M2; the sky location of the signal {9, (f>)\ the distance, 
D, to the signal; the coalescence time of the signal, to] 
the orientation of the binary, given by the inclination t, 
the polarization angle ?/; and the coalescence phase (po- 
We also assume that the sky location (6*, (f) of the sig- 
nal is known, thereby reducing the number of unkown 
parameters to seven. 

In the radiation frame, where the gravitational wave 
propagates in the e^-direction, the gravitational wave- 
form is given by 



h = /i+ea 



/ixe> 



(2.1) 



where 



ex = e^ 



(2.2) 



and the waveforms /i+,x depend upon seven parame- 
ters (All, M2,to, D, L^Tp, (po). The three remaining an- 
gles (t, ip, (po) give the relationship between the radiation 
frame and the source frame (in which ef lies in the direc- 
tion of the binary's angular momentum and ef along the 
separation between the binary components at to). Even 
for a known sky location, it is necessary to search a seven 
dimensional parameter space of signals. Naively covering 
this space with a grid of templates would be prohibitively 
costly 37 . However, the analysis is greatly simplified 
by the observation that the last four parameters enter 
only as amplitude parameters which can be analytically 
maximized over at minimal costj^ Specifically, the two 
polarizations of the waveform can be expressed as 



h+{t) = A^ho{t) + A^h^{t) 
h^{t) = A^h^{t) + A^h^(t) ^ 



(2.3) 



The two phases of the waveform are written as /iq and hjL . 
These depend upon the physical parameters of the system 
(in this case just the masses) as well as the coalescence 
time to^ A^ are constant amplitude terms and are given 



^ This was observed for the inspiral signal in I30| and indepen- 
dently for continuous wave signals in 1331 . 

^ This decomposition is actually valid for all binaries in which the 
plane of the orbit does not process. Thus, binary coalescence 
waveforms in which the spins are aligned with the orbital angular 
momentum can also be expressed in this form. However, for 
generic spin configurations, the orbit will precess and this simple 
decomposition is no longer applicable 
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explicitly as 



A" 
A^ 



cos 2(/)o cos 2?/; — ylx sin 20o sin 2-0 (2-4) 
A^ cos 2(/)o sin + Ax sin 2(^o cos 2-0 
— sin 200 cos 2-0 — A X cos 20o sin 2?/; 
— sin 20O sin 211; + Ax cos 20o cos 20; , 



where 



A, 
A, 



Do{l + cos^ i) 
~D 2 
Do 
D 



cos t , 



(2.5) 



and Do is a fiducial distance which is used to scale the 
amplitudes A^ and waveforms /io.f ■ Thus, the ampli- 
tudes A^ depend upon the distance to the source and the 
binary orientation as encoded in the three ang les U , tp, 
0o). For any set of values A^, the expressions (2.4) can 



be inverted to obtain the physical parameters, unique up 
to reflection symmetry of the system [36j . 

The gravitational waveform observed in a detector X 



+ polarization. The angle ip then describes the orienta- 
tion of the source with respect to this preferred radiation 
frame. 

Since we are considering CBC observed in ground- 
based detectors, the time that a potential signal would 
spend in the sensitivity band of any detector will be short 
(less than 60s for the initial detectors). Thus the change 
in the source's sky location over the observation time may 
be neglected, and the angles {0-^ , cj)^ , x'^ ) can be treated 
as constants. When working with a network of detec- 
tors, it is often useful to work in the geocentric frame. 
The location of the source {9, 0, %) is measured relative 
to this frame, and coalescence time is measured at the 
Earth's centre. In this case, the location and orientation 
of the detector X are specified by three angles, which we 
denote d'^ , and the detector response will depend upon 
six angles {a^ ,6,(t),x)- Then, the observed signal in a 
given detector 



h'^it) ^ F+(a^, 0, 0, x)/i+(t^) + Fx («^, 0, 0, x)h^ it"") , 

(2.11) 

where is the time of arrival of the signal at detector 
X, 



''J 



(2.6) 



where D^ denotes the detector response tensor. For an 
interferometric detector, the response tensor is given by 



D 



X 



X 



■.X\ 



(2.7) 



where the basis vectors and point in the directions 
of the arms of the detector. It is often convenient to re- 
express the gravitational wave signal observed in a given 
detector as 



h"" it) = F+ [9'' , 0-^ , )h+{t) + Fx (0-^ , 0-" , )/ix [t) , 

where the detector response to the two polarizations of 
the gravitational wave is encoded in the functions 



F+{9A.X) 
Fx{9,cp,x) 



1 



(1 + cos^ 0) cos 20 cos 2x 



cos sin 20 sin 2% , 

i (1 -I- cos^ 0) cos 20 sin 2% 

cos sin 20 cos 2x. 



(2.9) 



(2.10) 



These response functions depend upon the three angles 
{9-^ , (j)-^ , x'^) which relate the detector frame to the radi- 
ation frame: 0-^ and cj)'^ give the sky location relative to 
the detector, while x'^ is the polarization angle between 
the detector and the radiation frames. We have, some- 
what unconventionally, allowed for a polarization angle 
in transforming from source to radiation and radiation 
to detector coordinates. In what follows, we will often 
find it convenient to fix the angle x'^ by explicitly tying 
it to the detector (or geocentric) frame; for example, by 
maximizing the detector (or network) sensitivity to the 



^t - dt{d^ ,9,(t),x) 



(2.12) 



and dt gives the difference in arrival time of the signal 
between the geocenter and detector, for the given sky 
position. 

Combining the final expressions for the binary coales- 
cence waveform (2.3) and the detector response (2.11), 



we can express the gravitational waveform observed in a 
given detector as 

4 

/i^(i) = ^-4^(i^,V',0o,O/»^W (2.13) 



where the A^ are defined in (2.4) and h]^ are given by 



h^{t)=F^ho{t'') 
h^{t)^F^ho{t'') 
h^{t)^F^h.{t^) 
h^{t)^F^h^it^). 



(2.14) 



B. Multi detector binary coalescence search 

Matched filtering theory [39] provides a method for 
determining whether the signal /i(i,^), parametrized by 



^ We do not give the explicit formula for the response function 
dependent on the six angles x), as the expression is 

somewhat lengthy. It can be obtained by performing six succes- 
sive rotations to the detector response tensor to transform from 
the detector frame, via the geocentric frame, to the radiation 
frame. The calculation is detailed in |38| . 
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the time and other parameters ^, is present in a noisy 
data stream. The data output by a detector is 



(t) = n^(t) + /i-^(i,e) 



(2.15) 



where n^it) is the noise, taken to be Gaussian and sta- 
tionary. The noise (t) of the detectors is characterized 
by the noise power spectral density (PSD) (f) as 



("^(/)K(/')r)-^(/-/')^.(/)- 



(2.16) 



With this, we define the single detector inner product 
between two time series a(t) and b(t) as 



(a-^|&-^) =4Re 



oo ~X 



«-^( /)r(/r 



(2.17) 



Then, the likelihood ratio of there being a signal h present 
in the data is given by: 



Pis\h) 



and the log-likelihood can be written as 
\ogA^i.s\h)-^-ih\h). 



(2.18) 



(2.19) 



The likelihood ratio for multiple detectors is a straight- 
forward generalization of the single detector expression 
(2.18). Assuming that the noise in different detectors are 



independent, in the sense that 

{n''{f)[n^{fr)) - S^'^Sif ~ nS^if) , (2.20) 

the multi-detector inner product is simply given by the 
sum of the single detector contributions 



(a|b) 



X 



The multi-detector likelihood is given by 



1 



lnA=(s|h)--(h|h) 



(2.21) 



(2.22) 



This is the optimal statistic for signal detection in Gaus- 
sian noise: the larger the value of In A, the more likely 
that a signal is present. It is not, however, an optimal 
statistic if the noise is not Gaussian, as we will explore 
in more detail in sections IlIII and HVl 

Specializing to the case of binary coalescence, we can 
substitute the known waveform parametrization (2.13) 



into the general matched filter likelihood (2.22). The 
multi-detector likelihood becomes 



lnA==yl^(s|h^ 



1 



(2.23) 



The derivative of (2.23) with respect to A^ provides the 



values of A'^ which maximize the likelihood as 



(2.25) 



where, following [40], we take J\4^'^ to be the inverse of 
■M/ii/- We then define the coherent SNR via the maxi- 
mum likelihood as 



2 

Pcoh 



21nA|, 



(s|h^)X^''(s|h,) 



(2.26) 



It is not difficult to show that p^^j^ follows a dis- 
tribution with four degrees of freedom in the absence of 
a signal, and a non-central distribution (again with 
4 degrees of freedom) when a signal is present. See, for 
example, [ID] for more details. Furthermore, p^^j^ is now 
a function of only the waveform components and no 
longer the A'^ parameters. Thus four of the original seven 
waveform parameters have been analytically maximized, 
leaving three to be searched over. 

Calculating the maximized likelihood, as well as esti- 
mating the parameters A^ requires an inversion of the 
matrix A^^i/. CBC signals will spend a large number of 
cycles in the sensitive band of the detector and conse- 
quently the and ^ phases will be (close to) orthogonal. 
Since the frequency evolves slowly, the amplitudes of the 
two phases will be close to equalj^i.e. 



X\2 



Therefore, the matrix A4 simplifies to 



f A C 

C B 

A C 

\0 C B 



(2.27) 
(2.28) 



(2.29) 



where 



B = Y.x{<y''F^Y 

C = Y.x{'^''Fl^)i^^F^). 



(2.30) 



In this form the detection statistic is almost identical to 
the J-"-statistic for detecting rotating neutron stars, as 
described in [33], in the case where the neutron star has 
a small wobble angle. 



Dominant Polarization 

Since we have included a polarization angle in both the 
transformation between geocentric and radiation frame 
(x) and between radiation and source frame ('0), we have 



where the matrix AA^n is defined as 
M^,^ := (h^|h^) . 



(2.24) 



* Indeed, several CBC waveforms are generated directly in the 
frequency domain I41| . making these equalities exact. 
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the freedom to specify one of these without placing any 
physical restriction on the signal. The coherent SNR is 
further simplified by introducing a dominant polarization 
frame which renders Mfj,i, diagonal. 

Under a rotation of the radiation frame by an angle 
1 the detector response functions transform as 



F: 



X 



F 



DP.X 



= Ff cos 2x'''^ + F^ sin 2x"'^ (2.31) 
F^ ^ = -F^ sin2x°^ +F,^cos2x°^ 

The rotation through x^^ will have an identical effect on 
all detectors. Thus, there exists a polarization angle x^^ 
which satisfies 



C 



DP 



X 



^XpDP,X..XpDF^X. 



0. 



(2.32) 



This can be solved to give x as 



tan4x 



DP 



2E;.(^^^f)(^^^x^) 
Ex [{^''Ffr-ia^F^r] 



(2.33) 



This choice serves to diagonalize the matrix Ai. To 
uniquely determine x^^, we impose an additional re- 
quirement that the network be more sensitive to the + 
polarization than to the x polarization. The value of 
X^^ is a function of the detector network, the source lo- 
cation and waveform; in particular it depends upon F^ ^ 
and CT'^ . From now on, we assume that we are working 
in the dominant polarization frame and drop the DP su- 
perscript from our expressions. 

The concept of the dominant polarization frame has 
been introduced previously in un-modelled burst searches 
[421 - 144] . While the idea is very similar, the actual imple- 
mentation is somewhat different. 

In the case that both A and B are non-zero, i.e. that 
the detector has some sensitivity to both polarizations, 
the coherent SNR can be written, in the dominant polar- 
ization, as 



2 

Pcoh 



(2.34) 



(s|F+ho)V(s|F+h^)^ 

(F+ho|F+ho) 
(s|Fxho)V(s|Fxh|)' 
(FxholFxho) 

The coherent SNR can then be seen to arise as the 
quadrature sum of the power in the two phases of the 
waveform (0 and §) in the two gravitational wave polar- 
izations (4- and x). 



Network Degeneracy 

In many cases, a detector network is much more sensi- 
tive to one gravitational wave polarization than the other. 
In the extreme limit (e.g. co-located and co-aligned de- 
tectors such as those at the Hanford site) the network 
is entirely insensitive to the second polarization. In the 



dominant polarization frame, the network becomes de- 
generate as -B — )■ or equivalently 



X 



(2.35) 



Thus the network will only be degenerate if F:^ ~ for 
all detectors X. If the network is degenerate then it is 
easy to see that the detection statistic will be degenerate 
as well. In this case it is logical to remove the x terms 
from the detection statistic reducing it to 



2 

Pcoh 



(s|F+hoj 



(s|F+h; 



(F+ho|F+ho) 



(2.36) 



which is x^ distributed with two degrees of freedom. 

In this formalism the coherent SNR changes abruptly 
from (2.36) to (2.34). If there is any sensitivity, no mat- 
ter how small, to the x polarization, there is an entirely 
different detection statistic. This arises due to maximiza- 
tion over the A'^ parameters, allowing them to take any 
value. Thus, even though a network may have very little 
sensitivity to the x polarization, and consequently there 
be little chance of observing the waveform in the x polar- 
ization, this is not taken into account in the derivation. 
A possible modification is to place an astrophysical prior 
on the parameters (Z?, i, ip, 0o) and propagate this to the 
distribution of the A^ [IS] . This would provide a smooth 
transition from the degenerate to non-degenerate search. 



C. Comparison with Coincident Search 

The single detector search is a special case of the de- 
generate network (2.36) and can be written as 



P'x 



(2.37) 



A coincidence search requires a signal to be observed in 
two or more detectors, without requiring consistency of 
the measured waveform amplitudes in the different de- 
tectors. In many cases, coincidence searches have made 
use of different template banks in the different detectors 
[IKHIHT] and required coincidence between the recovered 
mass parameters [46] . A comparison with the coherent 
analysis discussed above is facilitated if we consider a co- 
incident search where an identical template is used in all 
detectors, as was done in an analysis of early LIGO data 
[47) . In this case, the multi-detector coincident SNR is 
given by 



2 2 

Pcoinc — 2^ PX 
X 



X 



X\2 



(2.38) 



This is not immediately comparable to the coherent SNR 
given in (2.34). However, both can be re-cast into similar 



forms by writing the coincident SNR as 



2 

r coinc 



EE 

X,Y i=0,f 



(^1^) [^""] (■' 



(2.39) 
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Similarly, the coherent SNR can be written as 



h. 



X.Y 1=0,1 

(2.40) 

where we have defined the orthogonal unit vectors (in 
detector space) , as 



X 



(2.41) 



The SNR of the coincident search (2.391 is simply the 



sum of all power consistent with the template waveform 
in each detector. The coherent SNR (2.40) makes use 



of the fact that gravitational waves have only two polar- 
izations to restrict the accumulated SNR to the physi- 
cal subspace spanned by and For a signal, the 
power will lie entirely in this subspace, while noise in the 
detectors will contribute to all components of the coinci- 
dent SNR. Thus, the coherent analysis obtains precisely 
the same signal SNR but reduces the noise background. 
Specifically, the coherent SNR acquires contributions 
from four noise degrees of freedom, while the coincident 
SNR has 2N noise degrees of freedom, where TV indicates 
the number of active detectors. For a non-degenerate two 
detector network, the coincident and coherent SNRs are 
equal as in this case /+"/+ + fx fx = ^'^^^ In the case 
where a network is sensitive to only one polarization, the 
coherent SNR is constructed solely from the f^ direc- 
tion and coherent SNR is distributed with 2 degrees 
of freedom. 

Finally, we note that restricting to the coherent SNR 
can help to separate transients from gravitational wave 
signals as those transients which do not contribute power 
to the signal space will be ignored. However, many noise 
transients will contribute to the coherent SNR and more 
active methods of removing them are required. These 
methods are the focus of sections Uni and IIVI 



D. Synthetic + and x detectors 

In the dominant polarization the coherent SNR is com- 
prised of separate + and x components, with no cross 
terms. We can go one step further and interpret the 
coherent SNR as arising from two synthetic detectors, 
one sensitive to only the -I- polarization and one sensitive 
to only the x polarization. These synthetic detectors 
are most easily formed by combining the "overwhitened" 
data streams o'^ from the various detectors, where 



S^if) 



(2.42) 



The overwhitened synthetic data streams are simply 



0+,x(/) = E^- 



and the power spectra for these overwhitened data 
streams are 



Using this, the un-whitened synthetic data streams are 
given a^ 

In terms of these synthetic detectors the detection statis- 
tic becomes 

2 _ is+\ho)l + is+\h^)l {sx\ho)l + {sx\h^)l 

Pcoh — ~ 



(^ol^o)> 



(2.46) 

where the subscripts x on the inner products denote 
the fact that the power spectrum of the synthetic detec- 
tors is used in their evaluation. 



III. SIGNAL CONSISTENCY BETWEEN 
DETECTORS 

As discussed in the introduction, due to the presence 
of non-Gaussian noise transients, it is essential to make 
use of signal consistency requirements within search al- 
gorithms to distinguish glitches from gravitational wave 
signals. 

Multi-detector analyses have made good use of signal 
consistency between detectors (see e.g. IMj). A partic- 
ularly powerful test is the use of a "null stream" [26] 
which, by construction, contains no gravitational wave 
signal. Many noise transients will contribute power to 
the null stream and can therefore be eliminated as candi- 
date events. In addition, requiring that the gravitational 
wave signal is recovered consistently between detectors 
can eliminate other noise transients, in our case this is 
equivalent to imposing restrictions on the recovered val- 
ues of the parameters, These two methods will be 
considered in turn. For matched filtering searches, re- 
quiring consistency between the observed signal and tem- 
plate waveform has also proven very powerful |24) . A full 
description of waveform consistency tests is presented in 
the next section. 



A. Null Stream Consistency 

The gravitational waveform consists of two polariza- 
tions. Thus for networks comprising three or more de- 
tectors it is possible to construct one or more null data 



X „x 



if)-. 



(2.43) 



^ There is some ambiguity in fixing the overall normalization of 
the synthetic detectors. We require that our synthetic detectors 
have the same sensitivity to the two polarizations as the original 
network did by requiring ^ ^ )+,x = T,xiP+,x'^^)'^ ■ 
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streams which contain no gravitational wave signal 
In the context of a coherent search for CBC signals, the 
null consistency tests arise quite naturally. In section 
|II C[ we noted that the coherent SNR can be thought 
of as a projection of the coincident multi-detector SNR 
onto a four dimensional signal subspace. The remaining 
dimensions in the coincident search do not contain any 
gravitational wave signal, but will be subject to both 
Gaussian and non-Gaussian noise. Thus, we can define 
the nuU SNR as 



and the un-whitened null stream is 



2 

Pcoinc 

EE 

X,Yi=0..^ 



2 

Pcoh 



(3.1) 



where 



XY 



6 



XY 



X fY 



(3.2) 



A gravitational wave signal matching the template h 
will provide no contribution to the null SNR, so we ex- 
pect that, for signals, this statistic will be distributed 
with (27V — 4) degrees of freedom. A noise transient that 
is incoherent across the data streams, may give a large co- 
herent SNR, but it is likely to also give a large null SNR. 
Thus requiring a small null SNR will prove effective at 
distinguishing incoherent noise transients from real grav- 
itational wave signals. Since the definition of the null 
SNR (3.1 1 makes use of the template waveform, gravita- 



tional waveforms which do not match the template h can 
contribute to the null SNR. 

We can go one step further and introduce synthetic 
null detectors in analogy with the synthetic -I- and x de- 
tectors. For concreteness, we describe the three detector 
case, but this can be extended in a straightforward man- 
ner to larger networks. To do so, we introduce the null 
direction 



^X ^ .XYZjYjZ^ 



(3.3) 



such that N-^^ = n-^nX . Then, the over-whitened syn- 
thetic null detector is 



X 



°^(/) = E^«''(/)- (3-4) 
The power spectrum of the null stream i^ 



X 



(n^)2 



(3.5) 



In this case, there is a normalization ambiguity. For the syn- 
thetic plus and cross streams, it was natural to require that the 
synthetic detectors have the same sensitivities as the original net- 
work. For the null stream this is not feasible as the network has 
zero sensitivity to a signal in the null stream, so we normalize 
such that {ho\ho)N = 1- 



j^XgX 



if) 



X 



■SnU). 



Finally, the null SNR can be written as 



N 



{ho\ho)N 



(3.6) 



(3.7) 



The null SNR described above differs from the multi- 
detector null stream formalism introduced in |26j and 
used by several other authors. A null stream is con- 
structed to be a data stream which contains no contri- 
bution from the and gravitational waveforms, re- 
gardless of the details of the waveform. To provide a 
concrete comparison between the null stream and null 
SNR, we restrict attention to a three detector network. 
The null stream is explicitly constructed as 



SNun{f)^exYzF^F^s^{f) 



(3.8) 



By comparing the null stream in (3.8) with the syn- 



thetic null detector (3.6 1, it is clear that these will gener- 
ically differ. To get an insight into the differences, con- 
sider a network with two co- located detectors A and B, 
with power spectra S^{f) and S^{f) respectively, and 
a third detector C which is sensitive to the other po- 
larization of gravitational waves. For this network, the 
null stream will be a combination of only the A and B 
detector data. The power spectrum of the null stream is 



while for the synthetic null detector, it is 



1 



1 



1 



(3.9) 



(3.10) 



Thus, if the power spectra of detectors A and B are 
identical, then the two null streams are also identical. In 
the extreme case that the sensitivity bands of the two de- 
tectors do not overlap at all, then there is no null stream 
(•SnuII — ^ oo). However, the null SNR need not vanish 
and is similar to a two bin version of the t^st de- 
scribed in section [IVD| Thus, it is possible to construct 
scenarios in which these two null stream formulations dif- 
fer significantly. 

For the most part, the power spectra of the ground 
based detectors are rather comparable. So, in general 
there will not be a significant difference between these 
two forms. There are advantages to both methods. The 
null stream is designed to cancel all gravitational wave 
signals from the data, thus making it more robust when 
the signal is not well known. However, by making use of 
the template signal, there are instances in which the null 
SNR provides a more powerful consistency test. Further- 
more, it has a computational benefit in that it does not 
require the production of a null stream — all manipu- 
lations are performed on the single detector SNR data 



8 



streams which are subsequently separated into coherent 
and null componenets. In practice we have found very 
little difference in performance, and choose to compute 
the null SNR (3.1) for computational simplicity. 



Finally, we note that both null stream formalisms will 
perform optimally only if the three detectors have simi- 
lar sensitivities. In the case where one detector is signifi- 
cantly less sensitive than the others, the null stream will 
generally tend to the data of that less sensitive detector. 
Also, the null formalisms described here will only com- 
pletely cancel a gravitational wave signal provided that 
the calibration of the data streams is accurate, any error 
in calibration will lead to a signal surviving in the null 
stream data. 
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B. Amplitude Consistency 

The four amplitude parameters A^, encoding the dis- 
tance to and orientation of the binary system, can take 
any values. Indeed, any set of corresponds to unique 
values of the distance, inclination angle, coalescence 
phase and polarization angle, up to symmetries of the 
system. However, some of these values will be signifi- 
cantly more likely to occur astrophysically than others. 
For example, the number of binary coalescence events 
is expected to be approximately proportional to star for- 
mation rate |11| and consequently should be roughly uni- 
form in volume. Thus, events are more likely to occur at 
a greater distance. Similarly, the gravitational wave am- 
plitude, at a fixed distance, is greater for face on signals 
than edge on ones, as is clear from (2.5). Thus, approxi- 



mately face on signals at a large distance are more likely 
signals than nearby, edge on ones. 

The distribution of noise events will follow its own 
characteristic distribution. For Gaussian noise, this dis- 
tribution can be evaluated and the signal and noise distri- 
butions incorporated into the ranking statistic IF . Non- 
stationarities in the data will again produce a different 
distribution of amplitudes. Specifically, the majority of 
transients are caused by a disturbance or glitch in a sin- 
gle detector with little or no signal in the other detectors. 
For networks with three or more detectors, this will typi- 
cally be inconsistent with a coherent signal across the net- 
work, leading to a large value of the null SNR. In certain 
scenarios, most notably for two detector networks, there 
will be a consistent set of values for the A^. However, 
these values carry the characteristic signature of a glitch. 
Specifically, the SNR contributions will typically be con- 
sistent with a nearby, close to edge on system {A^ ~ 0), 
with a very specific orientation to provide essentially no 
response in all but one detector. Thus, the glitch distri- 
bution of the A"" parameters, will be significantly differ- 
ent from the distribution expected for gravitational wave 
signals. In the remainder of this section, we explore the 
possibility of making use of the extracted A'' parameters 
to distinguish between glitches and signals. Unlike the 
null stream, amplitude consistency tests are available for 
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FIG. 1. The distribution of the recovered inclination an- 
gle plotted against coherent SNR for optimally oriented sig- 
nals (unfilled circles), uniformly distributed orientations (blue 
crosses) and simulated glitches (red circles). The top figure 
shows a network configuration where we are equally sensitive 
to both gravitational wave polarizations. The bottom figure 
shows a configuration where we are 5 times more sensitive to 
the -I- polarization than to the x . 



two detector networks. They should be especially useful 
in the case of the two 4km LIGO instruments, which have 
similar sensitivities to the majority of points on the sky. 

We have argued that the majority of gravitational wave 
signals will originate from (close to) face on binaries while 
the majority of noise transients will mimic (close to) edge 
on binaries. The recovered value of the inclination angle 
i should then serve to separate signals from noise. To in- 
vestigate this, we simulated a large number of simulated 
CBC signals and a large number of noise glitches; added 
Gaussian noise and plotted the recovered inclination an- 
gle in Figure [T] The glitches were generated as events 
with a large SNR in one detector coincident with Gaus- 
sian noise in a second detector. The signals were sepa- 
rated into two groups: the first with only face on binaries 
(I cost I = 1) and the second a uniform distribution over 
the two sphere (uniform in cos l and ip) of the binary ori- 
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entation. In both cases, they were distributed uniformly 
in volume and orbital phase. We also consider different 
network configurations, both containing two equally sen- 
sitive detectors. In the first case one detector is sensitive 
to + and the other to x polarization ; in the second case 
both detectors have strong and equal sensitivity to the 
+ polarization and weak but opposite sensitivity to the 
X polarization — rather typical for the Hanford, Liv- 
ingston network. For both sets of signals and choices of 
network, there is a clear distinction between signal and 
glitch distribution. However, there is a clear downwards 
bias on the recovered values of l. This can be understood 
by looking at the expressions for and Ax ■ For face on 
binaries, these will be equal but, in the presence of noise. 
Ax will be reconstructed to be somewhat smaller than 
A^. A relative difference of only 5% leads to a recov- 
ered inclination of 45° , so even for loud signals there can 
be large discrepancy between the actual and recovered 
inclination angle. 

Despite the difference in distribution between signal 
and noise, there is also a significant overlap of the popula- 
tions at low SNRs. Consequently, any threshold imposed 
on the recovered inclination angle is liable to either reject 
a fraction of signals or pass a fraction of glitches. It is, 
however, quite possible that knowledge of these expected 
distributions could be folded into the detection statistic 
in a Bayesian manner. 

We have found that using the observed SNR in the 
individual instruments to be a more effective discrimi- 
nator of signal and noise. To demonstrate the efficacy 
of such an approach, in Figure |2] we plot the single de- 
tector SNR as a function of the coherent SNR for the 
same population of glitches and the two classes of signals 
(face on and uniformly distributed orientation) described 
above. The glitches fall into two groups depending upon 
which detector suffered the glitch. Since our model de- 
tectors are equally sensitive, then on average one expects 
each detector to accrue 1 / \/2 of the coherent SNR. Even 
allowing for non-optimally oriented signals and the addi- 
tion of Gaussian noise, the signals follow this expectation. 
Only a small number of signals are found with SNRs in- 
consistent with the expected values, these are ones that 
have very specific orientations. Overall, the signal and 
glitch populations are very well separated, at least until 
the coherent SNR becomes rather small. 

The most effective strategy we have found is to require 
that all events have an SNR above 4 in the two most 
sensitive detectors in the network. The cut is illustrated 
in Figure |2] This strate gy removes the majority of glitch 
signals while having a negligible effect on the signal pop- 
ulation at large SNR. For lower SNR the signals which 
are lost due to this cut would be unlikely to be detec- 
tion candidates as Gaussian noise alone produces similar 
events. 
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FIG. 2. The distribution of single detector SNR plotted 
against coherent SNR for optimally oriented signals (unfilled 
circles), uniformly oriented signals (blue crosses) and glitches 
(red circles). The top figure shows a two equally sensitive de- 
tector network configuration where one detector sees only the 
4- polarization and the second detector sees only the x po- 
larization. The bottom figure shows a two detector network 
configuration where both detectors have strong and equal sen- 
sitivity to the + polarization and weak but opposite sensitiv- 
ity to the X polarization. This is meant to emulate a typical 
instance for the Hanford, Livingston network. The diagonal 
solid black line shows the expected SNR for the optimally ori- 
ented signals. The horizontal dashed black line indicates an 
SNR of 4. 

IV. COHERENT TESTS 

Data from gravitational wave detectors contain numer- 
ous non-stationarities due to both instrumental and en- 
vironmental causes. These non-stationarities, or glitches, 
typically do not match well with the CBC waveform. 
However, they often contain enough power that, even 
though the match with the template is poor, a large SNR 
is observed. In the previous section, we have seen how 
the use of various coherent consistency tests can mitigate 
this problem. Additionally, a number of other signal con- 
sistency tests have been implemented PH Hi] and used 
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in searches for CBC signals [H [571 HE]. These tests are 
all designed to eliminate glitches which have a different 
signal morphology than the template waveform. This is 
essentially done by testing whether the detector data or- 
thogonal to the signal is well described as Gaussian and 
stationary — for a glitch, there will be residual power 
which does not match the template waveform. These 
tests are commonly known as "x^ tests" as they con- 
struct a statistic which is distributed in the presence 
of Gaussian noise plus a signal matching the template 
waveform. If the data contains a glitch, the statis- 
tic will generally have a large value, thereby allowing for 
differentiation of signal from non-stationary noise. In 
this section, we briefly review the general formulation of 
X^ tests before presenting a detailed description of three 
such tests which have been implemented for the coherent 
search described in section |lll 



A. A general framework for x'^ tests 

Consider the data from a gravitational wave detector 
at a time t which has produced a large SNR when filtered 
against a template h{t). Generically, the data s{t) can 
be decomposed as 



s{t) = n{t) + Ah{t) + Bg{t) 



(4.1) 



where n{t) represents a Gaussian noise component, h{t) 
is the template waveform, g[t) is an additional non- 
Gaussian noise contribution to the data stream and A 
and B are amplitude factors. The glitch contribution 
g{t) is taken to be the power orthogonal to h{t) and both 
g{t) and hit) are normalized, so that 



{g\g)^l,{h\h)^l,{g\h)^0. 



(4.2) 



In order to construct a x^ test, we must introduce an 
additional set of waveforms T'. These waveforms are 
required to be orthonormal and orthogonal to h, 



(4.3) 



Furthermore, for the x^ test to be effective, the T* must 
have a good overlap with the glitch waveform g{t). 
The x^ discriminator is constructed as 



N 



(4.4) 



When the data comprises only signal plus Gaussian noise. 



i.e. B = in equation (4.1) 



N 



(4.5) 



and the statistic is the sum of squares of independent 
Gaussian variables with zero mean and unit variance. 



Thus the test is x^ distributed with N degrees of free- 
dom, with a mean and variance of 



{X^)^N, Var(x')-2iV. 



(4.6) 



This is true for any set of waveforms T* given the above 
assumptions. 

In the case where the data are not an exact match 
to the signal, we take both A and B non-zero, i.e. any 
signal or glitch can be decomposed into a part Ah{t) pro- 
portional to the template under consideration plus a sec- 
ond orthogonal contribution Bg{t). Clearly, for different 
glitches, the waveform g{t) as well as the amplitude fac- 
tor B will be different. In this case the x^ test takes the 
form 



N 
i=l 



[{T\nf + 2B{T\n){T\g) + B^{T^\gf^ 



This has a mean 



N 



{X^) = N + B^Y.^T\9f 



(4.7) 



(4.8) 



and a variance 



Var(x^) = 2Ar + 4B2 V(r|.g)2. (4.9) 



The x^ test is distributed as a non-central x^ distribution 
with N degrees of freedom and a non-centrality parame- 
ter El 



N 



x^b^Y^^tW 



(4.10) 



The challenge in constructing a x^ test is to select the 
basis waveforms such they have large overlaps with the 
observed glitches in the data. If this is done successfully, 
then any glitch producing a large SNR will also give a 
large value of x^, inconsistent with a signal in Gaussian 
noise. 

In many cases, there is some uncertainty in the tem- 
plate waveform. For example, the post-Newtonian (PN) 
expansion used in generating CBC waveforms is trun- 
cated at a finite (typically 3 or 3.5 PN [iS]) order and 
there will be differences between this analytically calcu- 
lated waveform and the one provided by nature. There 
are similar uncertainties in waveforms obtained from 
numerical relativity simulations 49 . Additionally, to 
search the full parameter space of coalescing binaries, 
a discrete template bank is used which allows for some 
mismatch between the templates and any potential signal 
within the parameter space [37]. Normally the template 
bank is created so that the mismatch is no larger than 
3% at any point in the parameter space. Finally, there 
are uncertainties in instrumental calibration [50j which 
will affect the match between signal and template. 
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We model these effects by parametrizing the signal as 
H{t) = A{^/l ~ h{t) + em{t)) , (4.11) 

where m{t) is the component of H that is orthogonal to 
h [{m\h) = 0] and e encodes the mismatch between signal 
and template in the sense that 



{H\h) 



(4.12) 



In most cases, it is reasonable to assume a mismatch 
of less than 5%. The obvious counter-example is when 
searching for highly spinning systems using non-spinning 
waveforms, see e.g. [5T1 [5^ . 



Since (4.11 ) is a special case of (4.1 ) it follows directly 
that the mean and variance of the test in the presence 
of a mis-matched signal are 



N 



ix^) =iV + A2e2^(r|m)^ 

N 

Var(x') = 2iV + iAh^ ^(r |m)'- (4-13) 



Since the SNR of the signal is proportional to A, the ex- 
pected value for a mis-matched signal increases with 
the strength of the signal. However, for mismatched sig- 
nals x^ oc e^A^ while for glitches oc and provided 
eA <^ B the two can be separated. See [23 for a more 
detailed discussion. 

When introducing the x? test, we assumed that the 
were orthonormal and orthogonal to the template wave- 
form h. In practice, this can be difficult to guarantee. 
The signal consistency tests discussed in the remainder 
of this section are constructed from gravitational wave- 
forms. If one picks a set of gravitational waveforms, i*, 
there is no guarantee that they will be either orthonormal 
or orthogonal to h. We can, at least, construct waveforms 
which are orthogonal to h by introducing 



f - {t'\h)h 



(4.14) 



While this ensures {h\T''') = it does not guarantee or- 
thonormality of the T\ (r*|TJ) = Thus this method 
will not produce a x^ distribution, and will instead form 
a generalized x^ distribution. The mean of the distribu- 
tion remains N but the variance is increased, 



Var(x^) = 27V + 2^{T\T^f 



(4.15) 



It has been found, however, that this does not present a 
significant obstacle to using these tests, especially as the 
thresholds are tuned empirically |53j . 



Multi-detector tests 

In section |llj we derived a coherent multi-detector 
search for coalescing binaries. The search involves fil- 
tering four waveform components against the multi- 
detector data stream. Our initial discussion of x^ tests 
was limited to the description of a single phase tem- 
plate waveform h and test waveforms T'. The extension 
to a two phase waveform has been described previously 
'24' and here we extend that to a four component wave- 
form across multiple detectors, as is appropriate for this 
search. We begin by noting that the four waveform com- 
ponents are orthogonal in the dominant polarization 
basis. They are, however, not generally normalized, as 



(h^|h,) = M^, = diag(A, B, A, B) , 



(4.16) 



where A and B are defined in (2.301. Thus, we first 
normalize so that 



(h^|h^) = 5f,u- 



(4.17) 



To construct a network x^ test, we require a set of 
(4-component) normalized, test waveforms i^. The com- 
ponents 



rpl 



1-E.fe|h.)2 



(4.18) 



constructed to be orthogonal to h^, are used in the x^ 
test. Thus, the coherent, multi-detector x? test is 



4 N 

/J, — 1 2 — 1 



(4.19) 



Provided the test waveforms are orthonormal, in the 
sense that 



(4.20) 



the distribution for a signal matching /i^ plus Gaussian 
noise will be x^ distributed with AN degrees of freedom. 
As for the single phase filter, we cannot always guar- 
antee (4.20) is satisfied, although it is relatively simple 



to ensure the four components of a given template are 
orthogonal. This means that the statistic will not, in 
general, be x^ distributed: The mean remains AN but 
the variance increases to 



N 



Var(x^)=87V+2^ ^ [{Tl\Tlf-5''5,.]. (4.21) 



i.j — l 1 



B. The coherent bank x^ test 

The bank x^ test was designed to test the consistency 
of the observed SNR across different templates in the 




FIG. 3. The distribution of the bank ^ test for a single 
template /i, with a bank of size 10. The plot shows the distri- 
bution of the bank veto calculated for every time sample in 
128s of simulated Gaussian data (with no signal present). In 
the case that the ten bank templates are orthogonal, the ex- 
pected distribution is with 40 degrees of freedom (shown as 
the solid black line). As can be seen, the actual distribution 
follows the expected one closely. 



FIG. 4. The single detector auto-correlation of a gravitational 
wave inspiral signal from a 1.4,1.4 solar mass binary neutron 
star. Both phases of the waveform are shown. 



C. The coherent autocorrelation x test 



bank at the time of a candidate signal. It was first de- 
scribed in |55] for the case of a single detector. A glitch 
will typically cause a high SNR in many templates across 
the bank, while a real signal will give a well prescribed 
distribution of SNR across the template bank. 

The bank makes use of other CBC templates as the 
waveforms to construct the x^ test. These iV templates 
are taken from different points across the mass space. 
In implementing the bank x^i we choose a fixed set of 
template waveforms V" which remain the same for every 
template h in the search template bank. The bank x^ 



statistic is then constructed following (4.18) and (4.191 



The test is most effective when the set of is close 
to orthogonal so we select templates which are well 
distributed across the mass space, ensuring the overlaps 
(T^jT^) are small for i ^ j. Figure |3] shows the dis- 
tribution of the bank for a single template filtered 
against Gaussian noise. The set of fixed bank waveforms 
consisted of ten waveforms distributed over the full mass 
parameter space. Using these waveforms, the deviation 
from a x^ distribution is negligible. 

For the bank x^ to be effective, glitches in the data 
must have a good overlap with a reasonable fraction of 
the templates t\ While, in general, it is difficult to pre- 
dict the composition of glitches in the data, it seems rea- 
sonable to assume that glitches which produce a large 
SNR for the template h will also have a good overlap 
with other waveforms in the template space. Thus, the 
set of templates which is spread across the parameter 
space is suitable. 



Filtering a gravitational wave template against data 
containing a matching signal produces a peak in the SNR 
at the time of the signal. Furthermore, there is a charac- 
teristic shape of this peak which depends upon the tem- 
plate waveform and also the noise power spectrum of the 
data. An example of this "autocorrelation" for a binary 
neutron star template is shown in Figure |4] A noise tran- 
sient in the data will produce a peak in the SNR but it 
will typically lack the characteristic shape produced by 
a genuine CBC signal. 

The "auto" x^ test was designed to test the consistency 
of the SNR peak [25] . It is a similar test to the bank x^ , 
but where the bank investigates consistency in SNR 
across the mass space, the auto tests for consistency 
of the SNR time series. The set of templates V" are chosen 
to be the original template h with time shifts SV" applied. 
The values of St^ are all unique and chosen to be of the 
same time-scale as the auto-correlation of the template 
waveform (typically 0.1s or less) and the duration of non- 
stationarities in the data, which is similar. 

In Figure [5) we show the distribution of the auto for 
a single template waveform filtered in Gaussian data. For 
this result, forty waveforms were used, equally spaced 
with a 1 ms spacing, and all with coalescence times prior 
to that of h. Thus, the auto x^ is testing the consistency 
of the SNR time series for 0.04 seconds prior to the SNR 
peak. The overlap (t'|t-'_) depends only upon the differ- 
ence St^ — 5P and Figure [l] shows clearly that a significant 
fraction of the overlaps are far from zero. Consequently, 
the auto x^ test has a distribution with a large deviation 
from a x^ distribution with AN degrees of freedom. 
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FIG. 5. The distribution of the auto test for a single tem- 
plate h, generated with 40 time shifted templates, with shifts 
between 0.001 and 0.04 seconds. The plot shows the distri- 
bution of the auto veto calculated for every time sample in 
128s of simulated Gaussian data (with no signal present). In 
the case that the forty time shifted templates are orthogonal, 
the expected distribution is with 160 degrees of freedom 
(shown in black). As can be seen, the actual distribution dif- 
fers significantly from this due to the non-orthogonality of the 
t'^ waveforms. 



D. The coherent x test 

The "standard" test originally proposed in [24^ has 
been used as a discriminator in many gravitational wave 
searches for CBCs. Given the template waveforms and 
the detector sensitivity, it is possible to predict the ac- 
cumulation of SNR as a function of frequency. By cal- 
culating the observed SNR contribution from a number 
of frequency bins, and comparing to the predicted value, 
one can construct a consistency test. 

Formally, given a template h which produced a can- 
didate signal with an SNR of p, calculate N non- 
overlapping frequency windows such that the expected 
SNR is p/N in each. Then, calculate the actual SNR 

in each of these frequency bins and compare with the 
expected value by calculating 



N 



(4.22) 



For a gravitational wave signal matching the template h 
plus Gaussian noise, this statistic will be distributed 
with A^ — 1 degrees of freedom. Written in the form (4.22) 



it appears different from the general case we discussed 
earlier. In it was shown that it can be re-expressed 



in the form (4.4). 

This test can be extended to coherent, multi- 
detector searches. Indeed, in [54], the construction was 
applied to a coherent search for continuous gravitational 
waves. Here, we present the extension to a coherent CBC 



FIG. 6. The distribution of the x'^ test for a single template h, 
split into 16 non-overlapping frequency bins. The plot shows 
the distribution of the test calculated for every time sample 
in 128s of simulated Gaussian data (with no signal present). 
The observed distribution of values shows (shown in grey) is 
an excellent match with the expected x^ distribution with 
sixty degrees of freedom (shown in black). 



search. First, define 



(4.23) 



to be the SNR contribution in the ith frequency bin to 
the SNRj^ The coherent statistic is then constructed 



as 



N 



(4.24) 



i — l fl—1 



As all the components are orthogonal it is easy to see 
that this statistic will be exactly distributed with 
AN — 4 degrees of freedom. One can interpret this as 
the sum of the single detector values for the ho and 
hiL waveforms in the synthetic -I- and x detectors. Figure 
[g] shows the distribution of the "standard" x^, using 16 
frequency bins. The distribution matches the expected 
with 60 degrees of freedom. 

An alternative approach to applying the t^st to a 
coherent search was proposed in |29| . This approach in- 
volves calculating the values for each of the active 
detectors and using these values to veto glitches. 



Strictly speaking the frequency bins for the F^. and com- 
ponents will be different because, as we have noted in equation 
the PSDs for the synthetic + and X detectors are not 
However, usually the difference between the two is small 
enough that it can be safely ignored to avoid computing twice the 
number of filters. Alternatively, in 1241 a method was presented 
for calculating the standard test using unequal frequency bins, 
that method could easily be incorporated into a coherent search. 



equal. 
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V. IMPLEMENTATION AND PERFORMANCE 
OF A COHERENT SEARCH 

In the previous sections we have presented a method 
for detecting gravitational waves from CBCs in a net- 
work of detectors. The coherent SNR described in |TT] 
is ideal for distinguishing signals in Gaussian data and, 
in sections |III| and |IV[ we have introduced a number of 
strategies for discriminating between signal and noise in 
non-Gaussian, non-stationary data. Here, we describe an 
implementation of the targeted, coherent search for grav- 
itational waves from CBC. In addition, we demonstrate 
its efficacy by performing test analyses of simulated data 
and real data taken from S4. 



we have made use of a bank generated with the initial 
LIGO design spectrum, with a maximum total mass of 
40Mq and a minimum component mass of IMq, these 
are the same values as used in . This method enables 
us to demonstrate the efficiency of this coherent search 
it is not the optimal solution. A simple improvement 
would involve placing a template bank appropriate for 
the (maximally sensitive) synthetic -I- detector defined 
in equation ( 2.44[ ). In many cases, a network is signif- 
icantly more sensitive to one GW polarization and, in 
these cases, this template bank would perform well. 



3. Coherent SNR and Null Streams 



Implementation of a coherent triggered search 
for CBCs 



Here, we describe the main steps by which the algo- 
rithms described in sections [ll HI and |IV| have been im- 
plemented. The analysis is available in the LIGO Scien- 
tific Collaboration Applications Library (LAL) suite [55] . 
and makes use of a large number of tools and methods 
previously implemented in that library. 



1. Analysis setup 

A targeted, coincident search for gravitational waves 
from CBC associated to GRBs has been implemented, 
and used in a search of S5 and VSRl data [H]. The co- 
herent search pipeline uses many of the same definitions, 
and much of the same architecture to determine the anal- 
ysis details. Specifically, "onsource" time is [-5,-f 1) sec- 
onds around the reported time of the GRB; this is when 
a gravitational wave signal would be expected [201 ES] • 
The noise background is estimated using 1,944 seconds 
of "offsource" data split into 324 trials of 6 second length 
each. These are used to calculated the significance of any 
event occuring in the onsource. As in Ref. |16j we impose 
a 48s buffer zone between the onsource and offsource re- 
gions. To obtain an accurate estimate of the detectors' 
power spectra, we only analyse data from a detector if it 
has taken at least 2190s of continuous data around the 
time of the GRB. Modulo this restriction, the coherent 
analysis is designed to make use of data from all detectors 
that were on at the time of the GRB. 



2. Template bank generation 

The problem of placing a template bank for a single de- 
tector has been extensively studied [571I57H55] . However, 
less thought has been given to the problem of placing an 
appropriate bank for a coherent analysis. Our current 
method is to use a template bank generated for one of the 
detectors in the network. In the results presented later. 



The data are first read in and conditioned using the 
methods and algorithms developed for the S4 search for 
post-merger ringdowns from CBCs [60l [61]. The data 
are downsampled to a frequency of 4096Hz and split into 
overlapping 256 second segments for analysis. The noise 
PSDs are calculated using the same method as in [6P] . 

Each template of the bank is filtered against the data 
from each detector to generate the single detector filters 
(s'^lft-^x) and sensitivities . The algorithms used are 
taken from the LAL FindChirp library ^B?, specifically 
those written to perform a search for spinning waveforms 
[52] using the physical template family (PTF) waveforms 



[63] _J The waveform templates are generated using the 
Taylor T4 post-Newtonian approximant [H]. The sin- 
gle detector filter outputs are shifted in time to account 
for the relative delays from the given GRB sky location. 
They are then combined to form the coherent and the 
null SNRs as described by ( |2.40| ) and A "trigger" 

is recorded at any time the coherent SNR is greater than 
6, and no louder event occurred in any template in the 
bank within 0.1 seconds. 



4- Calculating the tests 

The analysis calculates signal based vetoes in the same 
manner that it does the coherent SNR: The necessary 
single detector filters are constructed and then these are 
combined together to create the tests as described in 
the earlier sections. 

Calculating the "standard" test is computationally 
expensive. Therefore this veto is only calculated for a 
segment if there is at least one event within that segment 
with SNR larger than the threshold and with values of 
bank x^j auto and null SNR that do not immediately 
lead to it being dismissed as a glitch. 



* This choice stems from the desire to extend this search to incor- 
porate a single spin. This is particularly appropriate for NSBH 
binaries where the spin of the NS can be safely neglected. 
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FIG. 7. The distribution of SNR triggers in the offsource 
region plotted against time for an analysis of simulated Gaus- 
sian noise in the initial LIGO (HI, H2, LI) network. 



5. Simulated Signals 

The sensitivity of the analysis to gravitational wave 
signals from CBC is assessed by adding simulated signals 
to the data stream before performing the search. The 
analysis uses the architecture from [5D] to perform these 
simulations. A simulation is deemed to have been recov- 
ered by the analysis is there is an event within 0.1 sec- 
onds of the signal time; no attempt is made to guarantee 
a good match between simulated and recovered param- 
eters. While this does mean that, in principle, signals 
may be found due to a nearby glitch, the effect is mini- 
mal, particularly when considering detection candidates 
which are louder than all background. 

The results shown in the remainder of the section are 
based on a set of simulated signals comprised of binary 
neutron stars (component masses limited to be between 
1 and 3 Mq) at appropriate distances to be observable by 
the detector network. All of these binaries are oriented 
face on to the detectors. These choices are motivated by 
the fact that the current implementation is designed as 
a search for gravitational waves associated with GRBs. 

B. Analysis of simulated data 

The analysis was first run on simulated data for the 
initial LIGO network (HI, 112, LI). Data were simulated 
to be Gaussian and stationary, with no noise transients. 
Specifically, the coherent analysis pipeline was run on 
2190 sections of Gaussian data as if a GRB had occurred 
in the middle of the data stretch. This provides a bench- 
mark with which to compare results when running on 
real data. 

Figure [7] shows all the triggers produced by the pipeline 
in the off source time. The loudest event in the approx- 
imately 2000 seconds of off-source data has an SNR of 



10.0 
9.5 



9.0 




GPS time 



FIG. 8. The distribution of SNR triggers in the offsource 
region plotted against time for an analysis of a mock S4 GRB. 
The axes on the plot are chosen to be identical to those for 
figure [7| to make the plots easier to compare. The S4 data 
has a large number of non-Gaussian features. The largest of 
these peaks extends to a coherent SNR of 40, although non- 
Gaussian structure is visible at SNRs as low as 7. 

7.24. Ideally, the various signal consistency tests de- 
scribed previously will reduce the amplitude of the loud- 
est surviving event in real data to something similar to 
this. 



C. Analysis of real data 

A test analysis was performed on real data taken from 
S4. We chose an arbitrary block of 2190 seconds of data 
for which all three of the LIGO detectors were operat- 
ing and ran the analysis as if a GRB had occurred dur- 
ing this time. The simulated sky location of the GRB 
was (184. 623°, 42. 294°) in right ascension and declination 
respectively. For this chosen time and sky location the 
sensitivity of the HI and LI detectors were roughly equal 
and the H2 detector was half as sensitive as the other two. 



Coherent SNR 

Figure [8] shows the coherent SNR of triggers produced 
during the analysis of the S4 data. It clearly demon- 
strates that this data is not well characterized by Gaus- 
sian noise alone. A number of loud transients are present 
in the data which show up as short duration peaks of 
large SNR. The loudest of these has an SNR of almost 
40. If events were only ranked on SNR, a signal would 
have to be very loud to show against this non-Gaussian 
background. In addition to the loud peaks there are also 
a large number of quieter peaks that show up at all times 
in the data. If the search is to begin to approach the ef- 
ficiency of the Gaussian case then the signal consistency 



16 




Coherent SNR Coherent SNR 



FIG. 9. The distribution of the null SNR plotted against 
coherent SNR. The solid line at null SNR of 5.5 is the line 
above which triggers are vetoed. The dashed line at 3.5 is the 
line above which triggers are downweighted (see section V C I. 



tests must allow us to reject the quieter noise transients, 
as well as the loud ones. 

In the remainder of this subsection we demonstrate the 
performance of the signal consistency tests introduced in 
and 



sections III 



IV 



Finally, we formulate a detection 
statistic, which arises as a combination of the coherent 
SNR and signal consistency tests, and demonstrate that 
with it the majority of the non-Gaussian background can 
be removed. 



Null SNR 

Figure [9] shows the performance of the null stream for 
both simulated signals and background noise. The abil- 
ity of the null SNR to distinguish signal from noise is 
relatively poor in this example. The mock GRB analysis 
uses data from the two Hanford detectors and the detec- 
tor at Livingston. Thus, the null stream is derived from a 
combination of the HI and H2 detectors; the Livingston 
detector does not contribute. The loudest glitches dur- 
ing the time of this analysis originated in LI, and there- 
fore do not have significant null SNR. However, quieter 
glitches in the Hanford detectors at an SNR around 10 
do produce a large null SNR. Any trigger with a null 
SNR greater than 5 is eliminated from the analysis. 



Single detector SNR 

The most straightforward, and most effective ampli- 
tude consistency test we have found is the requirement 
of a single detector SNR greater than 4 in the two most 
sensitive detectors; In this analysis, the LI and HI detec- 
tors. Figure [TO] demonstrates that this is a particularly 
effective strategy for removing noise glitches. Triggers 
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FIG. 10. The distribution of single detector SNR for the more 
sensitive HI and LI detectors, plotted against coherent SNR. 
The top figure shows the HI SNR, the bottom figure shows 
the LI SNR. The horizontal line indicates SNR=4. Below 
this line triggers will be vetoed. The inclined dark gray line 
indicates the expected SNR of these face on simulated signals. 



arising due to glitches in the LI detector have large co- 
herent SNR but a negligible contribution from HI and 
are consequently discarded. 



X tests 



In section IV we introduced three tests designed to 
sep arate signals from noise glitches in the data. Figure 
|ll| shows the distribution of the bank for every time 
sample for a single template. This is directly comparable 
to Figure|3]which shows the same for Gaussian data. The 
deviation from the predicted distribution is due to the 
non-Gaussianity of the data. 

The distribution of bank and auto foi' both sim- 
ulated signals and noise triggers is shown in Figure 12 
Both of these tests are effective at separating the sim- 
ulated signals from noise transients. In order to quan- 
tify this, we make use of the newSNR formalism that 
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FIG. 11. The distribution of the bank test for a single 
template /i, with a bank of size 10. The plot shows the distri- 
bution of the bank veto calculated for every time sample in 
128s of data. The observed distribution is inconsistent with 
the expected result in Gaussian noise (the black curve). 



FIG. 13. The distribution of standard x test plotted against 
SNR. The results for the off-source data triggers are plotted 
in blue, with simulated signals in red. The dashed lines show 
contours of newSNR, defined in section |V C[ the solid line 
shows the line of newSNR = 6. Triggers with newSNR < 6 
are vetoed. 
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FIG. 12. The distribution of bank and auto test plotted 
against SNR. The results for the off-source data triggers are 
plotted in blue, with simulated signals in red. The solid line 
shows the line of newSNR = 6, defined in section [V C[ Triggers 
with newSNR < 6 are vetoed. 



is being used in the latest coincident searches for CBCs 
[64l [65] . The idea is to downweight the significance of 
noise triggers with large values relative to signals. 
This is achieved by introducing the "newSNR" : 



P, 
P 



X < "dof 



1 + ^ 

"do/ 



2 4/3 



/2 



(5.1) 



where ndof is the number of degrees of freedom of the 
test. For a signal, the mean value is one per degree 
of freedom and consequently the newSNR will be similar 
to the SNR. Noise transients with a large x^ value are 
significantly downweighted. The newSNR is calculated 
using both the auto and bank values. Any trigger with 
either an auto newSNR or bank newSNR less than 6 is 
discarded. The curves on Figure [12] show this newSNR 
threshold for the two tests. 

Finally, we turn to the standard test. As this is 
rather costly to compute, we only do so for triggers which 
have passed all of the previously described thresholds (on 
coherent, null and single detector SNR, the bank and 
auto x^). Figure 13 shows the distribution of the stan- 
dard x^ test for simulated signals and noise. The preced- 
ing tests have succeeded in removing the vast majority 
of non-Gaussian triggers from the data. A threshold of 
6 on newSNR serves to eliminate a few more. We have 
found that the standard x^ is the most effective at sep- 
arating signal from background so we also make use of 
it in the final ranking of events. Figure |13[ shows con- 
tours of constant newSNR which will be used in the final 
ranking. 
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Detection statistic 

In the preceeding discussion, we have imposed a num- 
ber of cuts on the initial candidate events produced by 
the analysis pipeline. Let us briefly recap those cuts: 

• Generate a trigger at any time for which p > 6. 
Only keep the loudest one in each 0.1 seconds. 

• Discard any triggers with > 5. 

• Discard any triggers for which pni < 4 or p^i < 4. 

• Discard any triggers for which Pnew < 6 for the 
bank or auto x^- 

Finally, we rank the remaining triggers based upon the 
newSNR calculated using the standard as well as the 
null SNR as: 

• Rank remaining triggers using a detection statistic 
Pdet given by 



Pdot 



Pncw 5 

Pncw 

PNuU — 2.5 



Pn < 3.5 

3.5 < Pn < 5.0 



(5.2) 
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FIG. 14. Efficiency of recovery of simulated signals in real 
data. The efficiency is shown for four different cases: i) signals 
found above an SNR of 7.24, the loudest background trigger 
in Gaussian noise (black solid line); ii) signals found above an 
SNR of 12.88, the loudest background trigger in real noise (red 
dashed line); iii) signals found with a value of the detection 
statistic ( |5.2[ ) above 7.41, the loudest background event; iv) 
signals found louder than all background in an Hl-Ll search 
(green dotted line). 



The length of time a CBC spends in the sensitive band 
of the detector varies greatly with the mass, and it has 
been found that the shorter, high mass templates are 
more susceptible to occuring with large SNR at the time 
of glitches . Also, the various signal consistency tests 
are less effective for these short templates. Therefore, we 
follow Ref. [55] and split the template bank into three 
regions based on the chirp mass of the template. The 
false alarm probability for a given trigger is based on a 
comparison of the detection statistic to off source triggers 
in the same mass bin. 



Performance of search 

The sensitivity of the analysis can be assessed by ex- 
amining the performance of recovering simulated signals. 
The signal population was BNS with optimal orienta- 
tion originating from the location of the fake GRB. For 
these purposes we will only consider an signal found if its 
associated trigger is louder than the loudest event that 
occured in the low mass bin during the offsource time 
— this ensures that the false alarm probability is less 
than one in 324 (the number of offsource trials). At this 
level, a candidate event starts to be interesting but real- 
isitically a false alarm probability closer to 10"'' or 10^^ 
would be required for a detection candidate. 

Figure [14] shows the efficiency with which simulated 
signals are recovered as a function of the distance. Re- 
call that the largest SNR recorded in the analysis of sim- 
ulated, Gaussian data was 7.24. We show the efhciency 
of the search at finding simulations with an SNR greater 



than this. If the data were Gaussian, or the signal con- 
sistency tests were able to remove all non-Gaussianities, 
then the real search would match this sensitivity. By us- 
ing only SNR to rank events, the sensitivity of the search 
is substantially worse than in Gaussian data. Although 
the offsource contained events with an SNR up to 39, the 
loudest low mass trigger had a coherent SNR of 12.88. 
This is almost double the loudest Gaussian noise and 
thus the sensitivity of the search is reduced by about a 
factor of two — illustrated by the 50% efficiency moving 
down from 20 Mpc to 10 Mpc. The signal consistency 
tests — null stream, tests and single detector SNR 
threshold — are designed to bring the sensitivity closer 
to the Gaussian case. The simulated signals are consid- 
ered found if they are recovered with a detection SNR 
greater than 7.41 (the loudest off-source event) and all 
signal consistency tests. The efficiency of the detection 
search is only about 10% less than it would be in Gaus- 
sian noise. The 10% loss in sensitivity can be attributed 
to having a slightly louder offsource event (7.41 rather 
than 7.24) and a small loss of efhciency due to the vari- 
ous signal consistency checks. 

To illustrate the benefits of a coherent search, we 
would like to compare the performance with a coinci- 
dence search. For the initial LIGO network, signals close 
to detection threshold would be unlikely to have passed 
any SNR threshold placed on the H2 detector. Specif- 
ically, for a signal with a coherent SNR of 7.5, the ex- 
pected SNR in H2 would be around 2.5 which is a pro- 
hibitively low threshold. Therefore it is reasonable to 
compare the coherent search to a two detector, Hl-Ll 
coincidence search. For a two detector search, the coher- 
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ent and coincident SNRs are equal, and the null stream 
test is not available. Consequently, the performance of 
two detector coherent and concident searches should be 
comparable. Therefore, we present the results of an Hl- 
Ll detector coherent search to give an indication of the 
performance of a coincidence search. The efhciency of 
this two detector search is also shown in Figure [14] The 
sensitivity of the two detector search is about a factor 
of 10% lower than the three detector coherent search. 
For a network of three approximately equally sensitive 
detectors, we would expect an even greater sensitivity 
improvement from employing coherent techniques. 

VI. DISCUSSION 

We have presented a formulation of a targeted coher- 
ent search for compact binary coalescences. For Gaus- 
sian noise, the coherent SNR would be ideal for separat- 
ing signals from the noise background. However, since 
data from gravitational wave interferometers is neither 
Gaussian nor stationary, we have also discussed a num- 
ber of methods of separating the non-stationary noise 
background from the signal population. These tests in- 
clude various tests, which were originally designed for 
use in single detectors. We have extended them to the 
network analysis and demonstrated their continued effi- 
cacy. Additionally, the coherent analysis allows for some 
additional tests which are not readily available in the 
coincidence case. The most significant of these is the 
null SNR which can be used to reject events which are 
not consistent with two gravitational wave polarizations. 
Additionally, we explored consistency tests between the 
recovered amplitudes of the gravitational wave and found 
that a simple SNR threshold on the two most sensitive 
detectors gave excellent results. 

The analysis described in this paper has been imple- 
mented and in the final section we showed results of a 
test run. This made use of the S4 data from the LIGO 
detectors. Although the data was far from Gaussian, af- 
ter the application of all of the signal consistency tests 
the results were remarkably close to what would be ex- 
pected in Gaussian noise. This analysis is available to be 
used in searches for GW inspiral signals associated with 
GRBs in more recent LIGO and Virgo data, such as S6 
and VSR2 and VSR3. 

There are a number of ways in which this analysis could 
be enhanced to broaden its use and increase its sensitiv- 
ity. First, a number of GRBs, particularly those observed 
by Fermi [57] and IPN jBS] are not localized sufficiently 
accurately that the error box can be treated as a point 
on the sky. Thus, it would be nice to extend this anal- 
ysis to allow for a region of the sky to be covered. This 
would require looping over the relevant sky points; incor- 
porating the correct detector sensitivities ^+,x and time 



delays. In principle, this would not greatly slow down the 
analysis as the majority of time is taken in performing 
the single detector filters and these would not need to be 
re-calculated. As well as looking at a patch on the sky, 
the analysis could be extended to cover the whole sky, 
as appropriate for an un-triggered search. This brings in 
a host of new complications which have been met and 
dealt with by other coherent search methods 44, 69^ . In 
order to obtain a good estimate of the background for an 
all sky, un-triggered search we would need to implement 
background estimation and time shifting the data would 
likely be the best way to do this. 

Since GRBs are thought to be rather tightly beamed, 
it is reasonable to take them as being face on, or close 
to. In this case, the gravitational waves are circularly 
polarized and there is, in effect, only a single polariza- 
tion. This opens the possibility of limiting the signal 
space to just this one polarization and adding an extra 
"null" test. Alternatively, it should be possible to per- 
form a Bayesian marginalization over the astrophysically 
expected distributions of the various parameters. 

The progenitors of short GRBs are thought to be BNS 
or NSBH. The search we have described is ideal for the 
BNS case as the spins of the neutron stars are unlikely 
to have a significant effect on the waveform. However, 
when one of the components of the binary is a black hole, 
the spin could be large. Furthermore, the mass ratio is 
likely to be relatively large. In this case, the spin of the 
black hole can have a significant effect on the observed 
waveform [63^. Consequently, we would like to extend 
this search to incorporate spin effects. The infrastruc- 
ture described in this paper can already accept spinning 
waveforms, but the implementation of signal based ve- 
toes proves somewhat more complex. Work is underway 
on this [70]. 
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